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We calculate the hadron spectra in high-energy pp, pA and AA collisions within a pQCD 
^ ' parton model. Analyses of experimental data show that the modification of the spectra 
^ ■ in pA and AA as compared to that in pp collisions is consistent with multiple parton 
Q . scattering scenario and is dictated by the transition from absorptive soft interaction at 
i low Pt to incoherent hard parton scattering at high pt- This analysis not only sheds 
new light on the limits to the physics analysis based on thermal fire-ball models but also 
0> ' provides a more quantitative baseline results on which one can study the effects of parton 
energy loss in the hadron spectra at high px- 
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1. Introduction 



Large-pr partons or jets are good probes of the dense matter formed in ultra-relativistic 
heavy- ion collisions [0] . The study of parton energy loss can shed light on the properties 
of the dense matter in the early stage of heavy-ion collisions. Large pt single-inclusive 
', particle spectra in nuclear collisions are sensitive p[ to parton energy loss. It is also a 
crucial test whether there is any thermalization going on in the initial stage of heavy-ion 
collisions. 

At low Pt the pQCQ parton model becomes invalid and other alternative approaches 
like thermal fire-ball models have to be used, from which one can extract the freeze- 
out temperature, collective radial fiow velocity and chemical potential 0. Apparently, 
these thermal fire-ball models cannot be applied to describe hadron spectra at large pt- 
Therefore, it is very important to investigate how well a pQCD parton model can describe 
hadron spectra in pp collisions and their modification in pA and AA collisions and where 
the transition happens between hard and soft hadron production. In particular, the 
impact-parameter or A dependence of the spectra may be unique to distinguish the parton 
model from other thermal fire-ball or hydrodynamic models. One can then at least make 
a quantitative conclusion about the validity of different models at different pt range. 
The values of temperature and flow velocity extracted from a flre-ball model analysis, for 
example, will have to be looked at with caution and knowledge of limitations. 



*This work was supported by the Director, Office of Energy Research, Division of Nuclear Physics of 
the Office of High Energy and Nuclear Physics of the U.S. Department of Energy under Contract Nos. 
DE-AC03-76SF00098. 
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2. Hadron Spectra in pp Collisions 

In a pQCD parton model, the inclusive particle production cross section in pp collisions 
is given by 



dyd?pT 



j dXadXb(fkaTd^kbTgp{kaT,Q^)gp{hT,Q^) 

abed 

fa/p{Xa, Q'^)fb/p{Xb, Q^) ^'''^ ^ ^(^^ (^) 

where D^^^{zc, Q"^) is the fragmentation function of parton c into hadron h as parame- 
terized in from e'^e~ data, Zc is the momentum fraction of a parton jet carried by a 
produced hadron. We choose the momentum scale as the transverse momentum of the 
produced parton jet Q = Pt/zc- We also use a factor K ^ 2 (unless otherwise specified) 
to account for higher order QCD corrections to the jet production cross section. 

One normally assumes the initial kx distribution gN{kT) to have a Gaussian form. In 
this study we relax the Gaussian form by assuming a variance which depends on px leading 
effectively to a non-Gaussian distribution, 

{kl)N{Q) = l(GeV') + 0.2a,(g2)g2_ ^2) 

The parameters are chosen to reproduce the experimental data at around SPS energies. 

Shown in Fig. |1| are our calculated spectra for charged pions as compared to the experi- 
mental data for p + p collisions at i?iab =200 GeV. Without the initial kx smearing the 
calculations significantly underestimate the experimental data. This is because the QCD 
spectra are very steep at low energies and even a small amount of initial kx could make a 
big increase to the final spectra. As the energy increases, the QCD spectra become flatter 
and small amount of initial kx does not change the spectra much. Such parton model 
calculations fit very well at energies from ^/s = 20 — 1800 GeV for pp and pp collisions 
[^]. The Pt dependence of -n^ j-n^ ratio is another indication of the dominance of valence 
quark scattering in large pt hadron spectra at SPS energies. 

3. pA and AA Collisions 

We assume that the inclusive differential cross section for large pt particle production is 
still given by hard parton-parton scattering, except that the initial transverse momentum 
/ct of the beam partons is broadened. Assuming that each scattering provide a kx kick 
which also has a Gaussian distribution, we can in effect just change the width of the initial 
kx distribution, 

(4)a(Q') = (4)iv(Q') + ^\Q^){yA{h) - 1). (3) 

The broadening is assumed to be proportional to the number of scattering v/S^ the 
projectile suffers inside the nucleus. We will use the following kx broadening per nucleon- 
nucleon collision. 
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Figure 1. Single-inclusive pion spectra in 
p + p collisions at i?iab = 200 GeV in pQCD 
parton model with (solid) and without (dot- 
dashed) intrinsic /ct- Data are from Ref. ||^. 



Figure 2. 7r° spectra m. A + B collisions 
at the CERN SPS energies in pQCD parton 
model with (solid) and without (dashed) kx 
broadening. Data are from Refs. [0,^. 



The pt dependence of the broadening reflects the fact that the distribution of soft kx kick 
for each scattering does not necessarily have a Gaussian form. 

The above calculation has been shown to reproduce the nuclear modification of the 
hadron spectra in pA data very well |P. Shown in Fig. ^ are the calculated inclusive 
spectra for produced tiq'ui A + B collisions. The pQCD parton model calculations with the 
kr broadening due to initial multiple scattering (solid lines) agree with the experimental 
data (WA80 and WA98) []^,^ well very at px above 2 GeV/c. No parton energy loss has 
been assumed in the calculations. The dashed lines are the spectra in pp collisions at 
the same energy multiplied by the averaged number of binary NN collisions given by the 
nuclear geometrical factor. The difference between the solid and dashed lines is simply 
caused by effects of kx broadening and nuclear modification of parton distributions inside 
nuclei. 



4. A Scaling of Hadron Spectra 

According to the pQCD parton model, the hadron spectra at large pt should scale with 
the number of binary nucleon-nucleon collisions if no nuclear effect is included. So if one 
defines a ratio, 

p / X _ da\B/dy£pT , . 

Rab[Pt) = TTT V , /r ,2 — 

{Num^ryjaa^p/dyd^PT 

between spectra in AB and A^A^ collisions normalized by the averaged number of binary 
collisions (Abinary), the ratio will be approximately one for spectra from hard parton 
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collisions. Because of absorptive processes, low pt particle production, which can be 
considered as coherent over the dimension of nuclear size, has much weaker A-dependence. 
In the wounded-nucleon model, soft particle production is proportional to the average 
number of wounded nucleons, the above ratio will become 

Rab ~ 0.5^/ A'/' + l/B'/'') (6) 

So the ratio as defined in Eq. (^) will be smaller than one at low px and larger than 
one at large pt- Such a general feature has been found to be almost universal in both 
pA and AB collisions. One interesting feature from this analysis is that the transition 
between soft coherent interaction to hard parton scattering happen roughly around Pt=1-5 
GeV. This is also the place where hadron spectra in pp collisions start to deviate from a 
pure exponential form. One can then expect that for spectra above this value of pt the 
underlying mechanism of hadron production is dominated by hard processes. 

At even higher px, the effect of multiple scattering becomes less important, so the ratio 
Rab will approach to 1 again (higher twist effect should be suppressed by l/p^), as shown 
by the pA data 0. However, at SPS energy, such a feature cannot be be fully revealed 
because of the kinetic limit. One will then only see the initial increase of the ratio due to 
the transition from soft to hard processes. Such a change of spectra from pp to pA and 
AA collisions in a limited kinetic range looks very similar to the effect of collective flow 
in a hydrodynamic model. However, models motivated by parton scattering have definite 
A-dependence of such a nuclear modification. Therefore one should take caution about 
the values of temperature and flow velocity extracted from such a fire-ball analysis of the 
spectra, especially if one has to reply on the shape of the spectra in the pt region around 
1 GeV. 

As one can also observe that the comparison of parton model calculation and the 
experiment data does not shown any evidence of parton energy loss 0. This implies that 
the life-time of dense partonic matter could be shorter than the mean free path of the 
propagating parton. It also shows that the dense hadronic matter which has existed for 
a period of time in the final stage of heavy-ion collisions does not cause any apparent 
parton energy loss or jet quenching. If one observes a dramatic suppression of high pt 
hadron spectra at the BNL RHIC energy as predicted |l|J^,|l^ , then it will clearly indicate 
an initial condition very different from what has been reached at the CERN SPS energy. 
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